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ABSTRACT 

An  analytical  model  of  a  jet  injected  normally  from  a 
flat  plate  into  a  uniform  crossing  flow  was  modified  to 
provide  an  improved  method  of  predicting  the  interference 
effects  arising  from  the  complex  flow  fields  induced  by  ship 
bow-thrusters.   This  model  is  an  extension  of  previous  work 
based  upon  a  description  of  the  jet  as  a  series  of  distrib- 
uted vortices.   The  analysis  takes  into  account  the  position 
of  the  effective  source  of  the  jet  and  the  blockage  due  to 
the  presence  of  the  jet  in  the  crossflow.   Improvement  in 
the  predicted  pressure  distribution  on  the  plate  was  obtain- 
ed by  adopting  a  new  formula  for  determining  the  effective 
jet  source  which  more  closely  approximates  the  experimental 
results.   Further  improvement  resulted  from  the  incorpora- 
tion of  a  symmetric  foil  shape  to  simulate  the  presence  of 
the  jet  and  its  highly  turbulent  wake.   The  accuracy  of  the 
model  was  evaluated  for  representative  jet-to-crossf low 
velocity  ratios.   Good  agreement  with  experimental  results 
was  achieved  for  large  portions  of  the  interaction  field 
with  particular  improvement  in  the  near-field  and  down- 
stream regions. 
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I.   INTRODUCTION 

Since  man  first  set  sail,  ships  and  boats  have  been 
steered  by  stern-mounted  rudders  or  rudder-like  devices. 
The  bow  has  been  left  relatively  uncontrolled.   As  the  size 
of  ships  increased,  they  became  more  difficult  to  maneuver 
in  close  quarters,  as  when  entering  port.   Tug  boats  were 
often  required  to  provide  lateral  thrust  at  both  ends  of  the 
ship.   This  was  a  satisfactory  solution  to  the  problem  for 
many  years.   It  did  not  provide  the  degree  of  independence 
desired  by  ship  operators,  however,  and  other  methods  were 
considered.   Although  experimental  lateral  thrust  units  were 
used  in  the  mid-nineteenth  century,  the  first  practical  unit 
was  not  installed  until  1955.   This  unit,  a  through- tunnel 
type  bow  thruster  installed  in  the  "Princess  of  Vancouver" 
proved  highly  successful  [1] .   Since  then,  bow  thrusters 
have  steadily  gained  in  popularity.   The  advantages  of  these 
units  include  providing  a  greater  degree  of  maneuverability, 
increased  safety,  increased  ability  to  operate  in  rough 
weather  at  oil  rigs  and  research  sites,  and  reduced  tug 
charges. 

In  the  past  twenty-five  years,  ports  have  become  increas- 
ingly crowded,  labor  and  fuel  costs  have  greatly  increased 
tug  boat  charges,  and  ships  have  been  built  that  are  much 
larger  and  less  maneuverable.   Thus  lateral  thrust  units 


have  not  only  become  economical,  but  they  are  standard 
equipment  for  many  new-construction  ships. 

As  more  bow  thrusters  came  into  use,  it  became  apparent 
that  although  they  work  well  at  zero  or  very  low  forward 
speed,  their  effectiveness  drops  considerably  as  the  ship's 
speed  is  increased  [2] .   In  view  of  the  large  number  of 
these  units  currently  installed,  there  is  suprisingly  little 
information  available  concerning  thruster/hull  interactions. 

Soon  after  the  first  unit  was  installed  in  1955,  it  was 
learned  that  the  operation  of  a  tunnel-type  thrust  unit 
causes  a  low  pressure  area  to  form  on  the  ship's  hull  [2] 
(see  Fig.  1) .   To  find  the  net  resultant  thrust,  this  suc- 
tion force  must  be  integrated  over  the  effected  area  and 
subtracted  from  the  thrust  of  the  unit.   When  a  thrust  unit 
is  operated,  water  surrounding  the  jet  is  entrained,  as  in 
an  eductor,  causing  a  low  pressure  area  to  be  formed  on  the 
ship's  hull  in  a  small  area  about  the  jet.   In  addition,  if 
a  ship  is  moving  forward,  a  plume  is  developed  aft  of  the 
thrust  unit  and  a  system  of  vortices  is  formed.   These  vor- 
tices induce  a  velocity  in  the  water  at  the  ship's  hull, 
causing  a  pressure  drop.   Additionally,  the  jet  itself  acts 
as  an  obstruction  to  the  flow  of  water  past  the  ship,  caus- 
ing the  flow  to  change  speed  as  it  passes,  thus  causing 
another  pressure  variation.   These  induced  velocity  and 
blockage  effects  are  so  severe  that  a  thrust  unit's 
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effectiveness  can  drop  as  much  as  fifty  percent  with  a  speed 
increase  of  only  three  knots  [1] . 

The  effects  of  bow  thruster  jet/hull  interactions  upon 
a  ship  motion  can  be  better  understood  by  considering  the 
equations  of  motion  for  a  ship  moving  in  a  horizontal  plane. 
The  general  forms  of  the  equations  are  [3] : 

m(u-vij/)  =  X    Surge  equation 

m(v+u^)  =  Y    Sway  equation 

I  $  =  N        Yaw  equation 

The  linearized  forms  of  these  equations  are  sufficient 
for  this  analysis  and,  making  the  ship's  center  of  gravity 
the  coordinate  system  reference  point,  they  become 
respectively: 


■X  (u-u  )  +  (m-x.  )u   =  X  (1) 

u    o        u 


■Y  v  +  (m-Y.)v  -  (Y  -mu  )r  -  Y-r  =  Y  (2) 

v        v       r   o      r 


-Nv-N.v-Nr+(I-N.)r=N  (3) 

v     v     r      z   r 


where: 


m  =  mass  of  the  ship 

u  =  velocity  in  the  x-direction  (initial  value  u  ) 

1  o 

v  =  velocity  in  the  y-direction 

•\)  =   yaw  angle 

r  =  angular  velocity  =  <ji 

I  =  mass  moment  of  inertia 

11 


X  =  excitation  forces  in  the  x-direction 

Y  =  excitation  forces  in  the  y-direction 

N  =  excitation  torque 

Dots  indicate  differentiation  with  respect  to  time. 

The  subscripts  u,  u,  v,  v,  r,  r,  denote  differenti- 
ation with  respect  to  the  subscript  variable. 

The  left  sides  of  these  equations  contain  the  velocity 
and  acceleration-dependent  terms  such  as  added  mass,  added 
inertia,  and  hydrodynamic  drag.   The  right  sides  implicitly 
represent  the  total  forces  in  the  x  and  y  directions  and  the 
total  moment  about  a  vertical  axis  through  the  center  of 
gravity.   Although  these  terms  are  used  to  represent  various 
excitation   forces  acting  on  a  ship  due  to  rudders ,  stabi- 
lizers, propellors,  etc.,  this  analysis  will  deal  solely 
with  the  forces  and  moments  caused  by  a  through-tunnel 
bow-mounted  thrust  unit. 

In  the  y-direction,  the  force  due  to  a  bow  thruster  may 
be  written  as: 


Y  =  Yfc  +  Y 
t    s 


where : 


Y   =  -m.u. ,  the  ideal  jet  thrust 

"^     J  J 


Y   =  the  induced  "suction" force  on  the  hull 


m.  =  the  mass  flow  rate  of  the  jet 


u.  =  the  jet  velocity  (positive  to  starboard) 

12 


Inspection  of  Fig.  2  reveals  that  -Y  /Y   is  nearly  linear 

in  the  I u  /u .   ratio  range  of  .1  to  .3.   Using  this  velocity 
3 

ratio  range  as  an  example,  an  approximation  of  the  right 
side  of  equations  (2)  and  (3)  can  be  made.   This  portion  of 
the  curves  is  of  interest  because  the  largest  decrease  in 
thruster  performance  takes  place  here.   Since  the  slope  and 
intercept  of  this  line  are  dependent  on  hull  shape,  an 
average  performance  curve  has  been  plotted  in  Fig.  3  to 
represent  a  "typical"  hull.   A  straight  line  having  a  slope 
of  5.3  and  an  intercept  of  -0.65  is  used  to  approximate  the 
curve  in  the  I u  /u .   ratio  (hereafter  R)  of  .1  to  .3. 


Thus: 


-Y  /Y   =  5.3R  -  .65 


Then  since 


Y.  =  -m.u. 
t      3    3 

Yc  =  m.u. (5.3R  -.65)  (4) 

and 

Y  =  -m.u.  (1.65  -5.3R)         .1  *  R  £  .3      (5) 

By  substituting  this  equation  into  equation  (2) ,  the 
motion  in  the  y-direction  is  explicity  coupled  to  the  motion 
in  the  x-direction  through  the  dependency  upon  R. 

The  total  turning  moment  N  is  the  moment  produced  by 
the  force  Y  acting  at  a  distance  x   from  the  center  of  gravity 
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This  characteristic  moment  arm  varies  non-linearly  as  the 
ship's  speed  changes  and  is  given  as: 


x^Y,.  +  x  Y 
t  t    s  s 

X   =  rz 

C  Y 


where: 


x   =  the  distance  from  the  center  of  gravity 

to  the  thrust  unit 
x   =  the  distance  from  the  center  of  gravity 

to  the  center  of  action  of  the  hull 

suction  force. 


Chislett  and  Bjorheden  [2]  found  that  x   varies  linearly 
with  ship's  speed.   The  curve  in  Fig.  4,  which  is  adapted 
from  their  data,  has  a  slope  of  -.7  and  an  intercept  of  1.0. 
Therefore: 


X   «  X. (1  -  .7R)  R  <  1.2  (6) 

s     t 

When  equations  (4)  and  (5)  are  substituted  for  Y  ,  Y   and  Y, 
the  location  of  the  center  of  action  of  the  combined  force 
system  is  given  by: 

x  =     (3'71R2  -  5-76R+  ^U      .1<R<.3    (7) 
C         (1.65  -  5.3R)       t 


and 


N  =  Yx 
c 
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Thus  the  yaw  equation  is  coupled  to  the  sway  and  surge 
equations.   Note  that  equation  (7)  is  only  valid  for  velocity 
ratios  between  .1  and  .3  for  this  particular  fictitious  hull. 

Although  this  analysis  has  been  greatly  simplified,  it 
demonstrates  the  complication  of  the  relationships  between 
turning  moments  and  ship's  speed  when  a  lateral  thrust  unit 
is  used.   It  can  be  seen  from  Fig.  2  that  model  tests  of 
each  ship  type  are  necessary  in  order  to  accurately  predict 
the  performance  characteristics  of  a  lateral  thrust  unit. 
With  a  thorough  understanding  of  the  fluid  mechanics  of  bow 
thruster  systems,  thrust  units  could  be  designed  to  reduce 
the  hull  suction  effect  and  thus  make  them  more  usable  under 
a  wider  range  of  conditions.   An  accurate  analytical  method 
for  predicting  jet/hull  interactions  would  be  useful  in 
supplying  this  information  without  lengthy  and  costly  model 
tests.   Unfortunately,  an  accurate  model  that  can  be  super- 
posed with  other  effects  is  not  available.   It  is  the  purpose 
of  this  thesis  to  propose  an  analytical  method  for  predicting 
the  jet/hull  interactions  of  a  through- tunnel  type  lateral 
thrust  unit. 
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II.   BACKGROUND 

As  demonstrated  in  the  preceeding  section,  the  ship's 
hull/thruster  interactions,  when  a  ship  has  way  on,  can 
seriously  degrade  a  bow-thruster ' s  effectiveness.   In  order 
for  naval  architects  to  design  a  bow  thruster  installation 
so  that  it  can  be  used  effectively  when  the  ship  is  moving, 
a  complete  knowledge  of  these  effects  is  required.   Unfortu- 
nately, only  a  few  waterborne  tests  have  been  performed. 
These  experiments  (Fig.  2) ,  while  clearly  demonstrating  the 
overall  effects  on  turning  ability,  do  not  provide  informa- 
tion on  the  actual  pressure  distribution  on  the  hull.   For 

an  incompressible  flow,  the  jet  path  and  the  induced  flow 

2      2 
are  dependent  mainly  on  the  momentum  flux  ratio  p  .  U_,  /p  U^ 

J  -i 

and  are  relatively  independent  of  the  Reynolds  number  [4]. 
Thus  an  air  jet  injected  into  a  crossing  air  flow  will  induce 
a  reaction  similar  to  that  of  the  water  jet  of  a  thrust  unit 
at  the  same  velocity  ratio.   There  have  been  a  number  of 
studies  of  air  jets  emitted  from  a  flat  plate  into  a  crossing 
air  flow  in  connection  with  the  design  of  VSTOL  (Vertical  and 
Short  Take  Off  and  Landing)  aircraft.   Notable  studies  have 
been  made  by  Bradbury  and  Woods  [4]  and  Fearn  and  Weston  [5]. 
While  the  results  of  these  experimental  studies  are  use- 
ful to  the  naval  architect  when  making  an  initial  design 
intended  to  improve  thruster  efficiency,  they  provide  no 

16 


method  for  evaluating  the  resultant  configuration.   Lengthy 
and  expensive  model  or  shipboard  tests  are  required  for  this 
purpose.   An  accurate  and  versatile  analytical  model  of  the 
flow  interaction  effects  would  provide  a  quick  and  inexpen- 
sive means  of  evaluating  new  designs,  thereby  reducing  the 
required  number  of  model  tests.   Some  empirical  and  semi- 
empirical  models  have  been  developed.   While  they  can  predict 
the  pressure  distributions  on  a  plate,  they  cannot  be  easily 
superposed  with  other  effects  such  as  a  secondary  jet. 
Because  a  fluid  jet  injected  into  a  crossing  stream  is  a 
turbulent,  three  dimensional,  highly  non-linear  problem, 
empiricism  to  some  degree  has  been  used  in  all  models  to 
date.   In  this  study  an  effort  has  been  made  to  eliminate  or 
at  least  reduce  the  empirical  inputs  to  existing  analytical 
models  while  improving  their  accuracy. 

A  review  of  available  models  revealed  that  Wooler's  formu- 
lation [6]  minimized  empiricism,  requiring  only  an  experi- 
mentally determined  jet  trajectory.   Hence  his  model,  as 
modified  by  Waterman  [7] ,  was  selected  as  a  basis  for  further 
investigation . 

In  order  to  visualize  the  interaction  field  and  determine 
the  accuracy  of  an  improved  model,  once  formulated,  accurate 
experimental  results  were  required.   A  review  of  available 
studies  revealed  that  the  pressure  contours*  for  a  given 


* 
In  this  thesis,  values  of  local  pressures,  P,  are  expres- 
sed in  the  form  of  pressure  coefficients,  C  ,  where: 

P-P  p 

C   =  **2 

p  Hp   tr 


velocity  ratio  vary  considerably  from  one  report  to  the  next, 
as  illustrated  in  Fig.  5  and  6.   This  difference  is  especially 
noticeable  in  the  downstream  region  for  values  of  |u./U 
(hereafter  VR)  of  about  4 ,  and  is  present  to  some  degree  at 
all  velocity  ratios.   Figure  5  seems  to  imply  that  the  induced 
flow  may  not  be  solely  a  function  of  velocity  ratio  since  the 
different  studies  were  not  made  at  the  same  freestream 
velocities.   However,  considering  studies  made  with  the  cross- 
flow  velocity  nearly  equal,  as  in  Fig.  6,  it  is  seen  that  the 
discrepancies  are  more  probably  due  to  variations  in  the 
equipment  and  test  procedures  used.   This  serves  to  point  out 
the  uncertainty  of  the  available  experimental  results.   The 
results  of  Fearn  and  Weston  [5]  agree  well  with  those  of 
other  researchers  over  a  broad  range  of  velocity  ratios. 
Additionally  their  report  was  the  most  complete  of  those 
reviewed.   It  contains  not  only  the  usual  plots  of  pressure 
contours,  but  also  raw  data,  time  variation  of  pressure 
information,  and  a  comparison  with  other  reports.   For  these 
reasons,  the  results  reported  in  [5]  were  selected  as  the 
basis  for  determining  the  accuracy  of  the  analytical  model. 
Before  attempting  to  develop  an  improved  model,  a  thor- 
ough review  of  Wooler's  and  Waterman's  work  was  made  in 
order  to  understand  their  approaches  and  verify  their  results. 
A  rederivation  of  these  basic  formulations  will  not  be  pre- 
sented here  as  it  can  be  found  in  references  [6]  and  [7] .   It 
was  noted  that  Waterman's  adoption  of  a  conforming  vortex 
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model  was  a  considerable  improvement  to  Wooler's  formulation. 
In  the  conforming  vortex  model  the  vortices  are  constrained 
to  lie  within  the  jet  rather  than  being  tangential  to  the  jet 
axis  as  in  Wooler's  model.   By  making  this  alteration,  the 
horseshoe  vortices  are  made  to  conform  to  the  curvature  of 
the  jet  described  by  the  empirical  trajectory  equation  (see 
Fig.  7) : 

|  =  B[cosh(Z/Bd)  -  1]  (8) 

Where : 

x  =  the  distance  down-stream  from  the  jet  center, 
d  =  the  jet  diameter. 

z  =  height  of  the  jet  axis  above  the  plate. 

2 
B  =  trajectory  parameter  =  . 19(VR)  . 

The  interference  velocity  at  an  arbitrary  point  on  the 
plate  is  determined  by  applying  the  Biot-Savart  Law.   The 
trailing  vortices  are  apDroximated  as  straight  line  segments. 
As  these  segments  become  small,  the  linear  approximation  can 
be  made  to  conform  to  the  curvature  of  the  jet.   The  total 
interference  component  at  a  particular  point  due  to  the 
entire  jet  results  from  the  summation  of  all  the  elemental 
contributions . 

A  comparison  of  the  results  of  both  formulations  reveals 
that  Wooler's  work  more  closely  agrees  with  the  experimental 
findings  of  Bradbury  and  Woods  [4]  than  does  Waterman's, 
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even  though  the  trailing  vortices  in  Wooler's  model  do  not 
follow  the  jet  trajectory  as  in  the  actual  case.   This  result 
is,  in  part,  due  to  the  fact  that  in  Wooler's  formulation  the 
perturbation  velocity  at  a  particular  point  on  the  plate  is 
caused  by  vortices  acting  at  different  distances  from  the 
plate.   In  Waterman's  model  the  interference  velocities 
result  from  the  trailing  horseshoe  vortex  contributions  act- 
ing at  a  single  radius.   Thus  the  interference  velocities 
predicted  by  Waterman's  conforming  vortex  model  are  generally 
higher  and  hence  produce  lower  coefficients  of  pressure  than 
those  of  Wooler.   In  addition  to  modifying  Wooler's  model, 
and  noting  studies  by  Keffer  and  Baines  [8]  and  Schmidt  [9] , 
Waterman  introduced  the  concept  of  an  effective  origin  of  the 
vortex  sheet  at  a  distance,  Z  ,  above  the  plate  given  by 
the  furmula: 


ze  =  -^vr 


Applying  this  expression  to  the  conforming  vortex  model, 
Waterman  produced  the  same,  if  not  better  far-field  agreement 
(Fig.  8) .   However,  the  model  lacked  accuracy  in  the  near- 
field. 

Waterman  further  noted  that  the  discrepancies  observed  in 
the  upstream  region  can  be  at  least  partially  attributed  to 
the  blockage  of  the  freestream  flow  by  the  jet.   It  was  noted 
that  the  downstream  area  is  dominated  by  separation,  vortex 
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shedding,  and  turbulent  entrainment.   Waterman  interpreted 
the  blockage  effect  as  that  of  a  circular  cylinder  in  an 
irrotational  flow.   Using  potential  flow  theory  to  predict 
the  velocity  perturbations  in  the  uniform  flow  caused  by  a 
doublet  centered  at  the  jet  origin,  Waterman,  after  some 
numerical  experimentation,  arrived  at  an  estimate  for  an 
effective  jet  diameter  given  by  .96  V  vA  .   The  results  of 
this  model  (Fig.  9)  are  significantly  better  than  Wooler's 
model  in  the  upstream  region  and  far  field.   This  improve- 
ment was  accomplished  at  the  expense  of  introducing  two 
additional  empirical  relationships,  the  position  of  the 
effective  source  and  the  effective  jet  diameter. 
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III.   IMPROVEMENT  IN  THE  EXISTING  MODEL 

Evaluation  of  Waterman's  "fully-corrected  conforming  vor- 
tex model"  (Figs.  9  and  10)  revealed  that  the  correlation  of 
the  analytical  model  with  the  experimental  data  in  the  area 
just  forward  of  the  jet  is  adequate  at  each  of  the  velocity 
ratios  reported.   However,  his  model  is  accurate  in  the  far- 
field  only  for  the  velocity  ratio  of  8.0  and  in  the  near- 
field  to  the  sides  of  the  jet  for  a  velocity  ratio  of  4.0. 
The  correlation  is  poor  in  the  area  aft  of  the  jet  for  all 
velocity  ratios.   It  was  apparent  that  the  near  and  mid- 
field  regions,  especially  aft  of  the  jet,  required  the 
greatest  attention. 

In  order  to  gain  a  better  understanding  of  the  effects  of 
the  empirical  inputs  to  Waterman's  model,  it  was  run  a  num- 
ber of  times  using  various  values  of  the  constants  in  the 
trajectory,  effective  source,  and  effective  diameter 
equations.   It  was  learned  that  small  changes  in  the  trajec- 
tory parameter,  B  in  equation  (8) ,  cause  little  change  in 
the  overall  pressure  distribution.   The  correlation  with 
experimental  data  was  improved  at  some  velocity  ratios  while 

it  diverged  at  others  for  a  particular  value  of  B.   A  value 

2 
of  B  =  .19(VR)  ,  as  before,  resulted  in  the  best  overall 

correlation  of  the  pressure  distribution  on  the  plate. 

Changing  the  effective  diameter  of  the  jet  causes 
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major  changes  in  the  near-field  pressure  distribution,  but 
has  little  effect  in  the  far-field.   This  was  expected  con- 
sidering the  short  range  of  the  disturbance  caused  by  a 
doublet  in  a  uniform  flow  using  inviscid  relationships.   The 
model  proved  to  be  particularly  sensitive  in  both  the  near 
and  far-field  regions  to  changes  in  the  position  of  the 
effective  source.   Since  the  quality  of  the  model  depends 
so  greatly  on  the  effective  source  correction,  the  equation 
for  it  was  investigated  more  closely. 

Equation  (9) ,  the  empirical  formula  for  the  effective 
source  height,  was  formulated  by  Schmidt  [9]  who  refers  to 
Jordinson  [10]  and  Keffer  and  Baines  [8]  for  the  experimental 
data  used  as  a  basis  for  his  work.   The  Jordinson  paper  could 
not  be  obtained  for  review.   The  data  presented  by  Keffer 
and  Baines,  however,  does  not  fit  the  equation  presented  by 
Schmidt.   The  equation,  Z  /d  =  . 8  V  VR  fits  the  Keffer  and 
Baines  data  much  better  as  seen  in  Table  I. 

TABLE  I 

COMPARISON  OF  VALUES  OF  Z  /d 

e7 


VR       Schmidt 


4  .53 

6  .97 

8  1.51 


Experimental 

Keffer  & 

Baines 

.8  VW 

1.59 

1.60 

1.95 

1.93 

2.30 

2.26 
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However,  Z  /d  can  vary  from  . 54  V  VR  to  . 8  V  VR  and  still  be 
in  the  range  of  the  heights  given  in  Table  I.   Using  an 


effective  source  location  of  Z  /d  =  .8  V  VR  in  the  Waterman 
model  significantly  improves  the  correlation.   But,  if 
Z  /d  =  . 6  V  VR  is  used,  the  correlation  is  much  better  even 
without  the  effective  diameter  correction.   Thus,  by  estab- 
lishing the  effective  source  at  . 6  V  VR  diameters  above  the 
plate,  the  full  jet  diameter  can  be  used  in  the  model  and 
the  diameter  correction  is  eliminated.   See  Fig.  11  for  an 
example  of  the  improved  correlation  resulting  from  this 
change . 

Concentrating  on  the  area  aft  of  the  jet,  it  can  be  seen 
by  comparing  Fig.  8  with  Fig.  10  that  this  area  is  primarily 
affected  by  the  blockage  near  the  plate  due  to  the  jet.   In 
making  multiple  runs  of  Waterman's  model  and  varying  the 
parameters  of  the  vortices  within  the  plume,  it  was  confirmed 
that  they  have  little  effect  on  the  near-field  pressures. 
The  far-field  to  either  side  of  the  jet  is  dominated  by  the 
vortices  within  the  plume.   Therefore,  to  improve  the  model's 
accuracy  in  the  near-field,  the  analysis  centered  on  the  jet 
blockage  effect.   Wu  and  Wright  [11]  proposed  that  a 
fictitious  afterbody  exists  downstream  of  the  jet.   This 
afterbody  encompasses  the  highly  turbulent,  three-dimensional 
flow  caused  by  vortex  shedding  and  entrainment  of  fluid 
within  and  downstream  of  the  jet.   Outside  of  this  area,  the 
fluid  streamlines  are  fairly  smooth  and  close  to  those  of  an 

inviscid  flow.   With  this  concept  in  mind,  consideration  was 
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given  to  using  an  elongated  shape  such  as  an  ellipse  or  a 
symmetrical  foil  to  approximate  the  blockage  effect. 

In  order  to  establish  a  probable  shape  for  the  afterbody, 
the  coefficients  of  pressure  along  the  corssflow  axis  from 
the  forward  edge  of  the  jet  to  seven  diameters  aft  of  the 
jet  were  plotted  using  the  data  of  [5]  (Fig.  12) .   These 
pressure  contours  resemble  those  in  reference  [12]  and  [13] 
for  a  symmetrical  foil  in  a  uniform  flow.   Thus  a  symmetrical 
foil  was  selected  to  approximate  the  jet  and  the  afterbody 
near  the  plate. 

To  determine  a  length  for  the  foil,  it  was  assumed  that 
the  foil  should  end  where  the  coefficient  of  pressure  reaches 
an  undisturbed  (zero)  value.   It  was  not  expected  that  the 
coefficient  of  pressure  should  necessarily  rise  above  zero 
aft  of  the  jet  due  to  severe  dissipation  in  the  near-wake 
region  and  the  interference  effects  of  the  vortices  within 
the  plume.   Indeed  the  coefficient  of  pressure  leveled-off 
at  approximately  -0.1  for  all  available  data.   Therefore, 
-0.1  was  selected  as  a  first  approximation  of  the  value  of 
the  coefficient  of  pressure  used  to  determine  the  length  of 
the  foil.   The  corresponding  distances  from  the  jet  center 
to  the  trailing  edge  were  determined  to  be  6,2,  3.0,  and  2.5 
jet  diameters  for  the  velocity  ratios  of  3.9,  8.Q,  and  10.0, 
respectively. 

For  simplicity,  a  standard  Joukowsky  foil  was  selected 
since  it  can  be  easily  approximated  through  mapping  methods. 
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Referring  to  Fig.  13,  the  Joukowsky  transformation, 


Z  =  c  +  9/ 


is  used  to  transform  the  foil  in  the  z-plane  to  a  cylinder 
in  the  c-plane  so  that  the  foil  can  be  treated  as  a  doublet 
in  a  uniform  stream.   The  perturbation  velocities  due  to  the 
presence  of  the  foil  are  then: 


and 


u 
U 


=  Jt 


m 


U   +  a) 
(?  +  I) 
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-  1 
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u 
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Where 


u  =  the  surface  perturbation  velocity  in  the  direction 
of  the  cross-flow  due  to  the  presence  of  the  jet 
and  afterbody  as  modeled  by  the  foil. 

v  =  the  surface  perturbation  velocity  perpendicular  to 
the  cross-flow  due  to  the  presence  of  the  jet  and 
afterbody  as  modeled  by  the  foil. 

I  =  1/4  the  distance  from  the  forward  focus  of  the 
foil  (corresponding  to  the  point  (-1,0)  in  the 
C-plane)  to  the  trailing  edge. 

a  =  the  radius  of  the  transformed  circle.   This  term 
is  a  function  of  r,  the  distance  from  the  forward 
focus  of  the  foil  to  the  leading  edge. 
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Appendix  A  contains  a  complete  derivation  of  the  above 
transformation . 

The  thickness  of  the  foil  is  a  function  of  both  r  and  I 
As  described  above,  the  foil  length  is  a  function  of  the 
velocity  ratio  so  that  the  length  I    is  given  by: 


I  -  -L  do 

d    VR  { 


Initially  it  was  assumed  that  the  forward  focus  of  the 
foil  should  be  at  the  center  of  the  jet.   This  makes  r  equal 
to  one-half  the  jet  diameter.   When  substituted  for  the 
doublet  in  the  fully  corrected  model,  this  method  provided 
significant  improvement  in  the  downstream  area.   The  up- 
stream correlation  was  poor,  however,  since  the  foil  was  too 
thick  and  too  blunt  at  the  leading  edge.   It  was  necessary 
to  adjust  the  value  of  r  to  provide  a  thinner  and  more 
streamlined  body.   Next  a  foil  with  its  thickness  equal  to 

one  jet  diameter  was  tested.   Using  the  formula,  maximum 

2 

thickness/d  =  5.2(m/2.  -  .  7(m/2,)  )  [13],  where  m  is  a  function 

of  r.   It  was  found  that  if  r/d  =  .126./ 1,    where  I    is  given  by 
equation  (10),  the  thickness  is  nearly  equal  to  the  jet 
diameter  for  the  entire  range  of  velocity  ratios.   The  foil 
is  then  moved  forward  so  that  its  leading  edge  coincides 
with  that  of  the  jet.   This  configuration  (App.  B)  produced 
the  results  shown  in  Figs.  14,  15,  and  16  in  comparison  with 
experimental  values  [5]  for  velocity  ratios  of  4.0,  8.0,  and 
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10.0  respectively.   Radial  and  circumferential  pressure 
distributions  for  velocity  ratios  of  3.9  and  8.0  are  plotted 
against  the  experimental  findings  of  Fearn  and  Weston  [5] 
in  Figs.  C-l  through  C-52. 

Attempts  were  made  to  use  other  fictitious  shapes  such 
as  ellipses  and  ovals  as  well  as  distributed  source-sink 
combinations.   The  foil  shape,  however,  produced  the  best 
results.   So  the  simple  relationships  for  the  near-field 
blockage  effect  given  above  were  retained  in  the  final 
model. 

Reviewing  the  results,  it  is  seen  that  for  a  velocity 
ratio  of  8.0,  the  correlation  between  this  model  and  the 
experimental  data  is  very  close  except  in  the  area  of 
highly  turbulent  flow  immediately  aft  (within  two  diameters) 
of  the  jet.   This  is  a  region  of  uncertainty  for  all  velocity 
ratios  due  to  the  three-dimensional  flow  caused  by  entrain- 
ment ,  vortex  shedding,  and  other  highly-viscous  wake 
phenomena.   Since  this  area  is  small,  it  contributes  little 
to  the  overall  thruster/hull  interaction.   For  a  velocity 
ratio  of  3.9,  the  correlation  is  excellent  except  for  the 
mid-  and  far-field  regions  aft  of  the  jet  centerline. 

This  was  expected  since  it  was  determined  earlier  that 
the  blockage  effect  of  the  jet  has  little  effect  in  the 
far-field.   It  was  observed  that  if  r  is  such  that  it  causes 
the  foil  to  have  a  circular  shape  at  the  leading  edge  the 
correlation  in  the  area  forward  of  the  jet  is  improved. 

28 


However,  the  agreement  deteriorated  in  the  areas  to  the 
sides  and  aft  of  the  jet.   By  varying  r/d  as  a  function  of 
the  location  of  the  point  at  which  the  velocity  perturbations 
are  being  calculated,  the  shape  of  the  foil  can  be  altered 
so  that  it  is  nearly  circular  at  the  leading  edge  and  becomes 
thin  aft  of  the  jet  centerline.   This  configuration  worked 
better,  but  it  failed  to  significantly  improve  the  overall 
pressure  distribution  resulting  from  r/d  =  . 12d/£.   This 
observation  did  lead  to  the  idea  of  using  a  half-cylinder 
shape  for  the  jet  blockage.   That  is,  the  fluid  forward  of 
jet  "sees"  a  cylinder  but  once  past  the  jet  centerline 
the  influence  of  the  cylinder  is  rapidly  diminished.   The 
results  of  this  configuration  are  presented  in  Fig.  17  from 
which  it  can  be  seen  that  while  this  is  a  significant 
improvement  over  Waterman's  model,  it  is  not  as  accurate 
as  the  full  foil-shape  model.   It  is  concluded  that  further 
improvement  of  the  model  must  involve  a  more  sophisticated 
treatment  of  the  vortices  within  the  plume  as  well  as  the 
blockage  effect  of  the  jet. 


29 


IV.   CONCLUSIONS 

As  evidenced  by  Figs.  14,  15,  16,  and  C-l  through  C-52 
this  analysis  adequately  predicts  the  pressure  distribution 
on  a  flat  plate  from  which  a  jet  is  emitted  into  a  crossflow. 
Application  of  the  model  is  limited,  however,  to  velocity 
ratios  VR  =  3.9  or  greater.   This  is  due  to  the  method  by 
which  the  length  of  the  fictitious  afterbody  is  determined. 
Above  velocity  ratios  of  about  4,  this  length  decreases  for 
increasing  VR,  but  plots  of  the  coefficients  of  pressure 
along  the  crossflow  axis  similar  to  Fig.  12  indicate  that 
the  afterbody  length  should  decrease  also  for  values  of  VR 
less  than  about  4.   This  is  consistent  with  the  observations 
of  Fearn  and  Weston.   They  reported  that  an  anomaly  in  the 
pressure  distribution  in  the  wake  occurs  at  a  velocity  ratio 
of  3.9.   The  implication  is  that  the  rate  of  pressure  recovery 
in  the  wake  region  is  a  minimum  for  a  velocity  ratio  of  about 
3.9  [5], 

No  attempt  has  been  made  to  accurately  predict  the  pres- 
sure distribution  in  the  highly  turbulent  wake  of  the  jet. 
This  area,  having  a  width  of  about  one  jet  diameter,  extends 
downstream  from  the  jet  along  the  crossflow  axis.   The  pre- 
dicted values  of  the  plate  pressures  reported  here  cannot  be 
relied  upon  since  most  of  this  area  is  contained  within  the 
fictitious  afterbody. 
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Figure  18  shows  the  results  of  this  analysis  plotted  as 
the  ratio  of  the  induced  side  force  to  the  jet  thrust, 
-Y  /Y. ,  versus  the  velocity  ratio  VR.   These  results  are 
compared  with  NPL  pressure  measurements  on  a  flat  plate  and 
the  model  tests  of  English  [1]  and  Chislett  and  Bjorheden  [2] . 
The  slope  of  the  predicted  -Y  /Y   curve  is  approximately  5.5, 
whereas  that  of  the  experimental  data  is  about  5.3  for  both 
the  flat  plate  and  the  model  tests.   While  this  analysis 
predicts  somewhat  larger  values  for  the  overall  (integrated) 
interference  effects,  the  small  difference  between  these 
slopes  indicates  that  the  nature  of  the  relationship  between 
the  interference  effects  and  the  velocity  ratio  predicted  by 
the  analytical  model  is  correct.   Fearn  and  Weston  did  not 
report  overall  interference  effects  so  the  correlation  between 
their  work  and  the  NPL  experiments  is  not  known.   Since  the 
assumptions  concerning  the  shape  and  size  of  the  afterbody 
were  based  on  Fearn  and  Weston's  data,  the  degree  to  which 
this  analysis  is  in  error  cannot  be  determined  from  Fig.  18 
alone.   Even  with  these  deficiencies,  it  is  felt  that  this 
analytical  model  adequately  predicts  the  pressure  distribu- 
tion on  a  flat  plate  and  can  be  used  as  a  first  approximation 
for  design  purposes. 
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V.   RECOMMENDATIONS  FOR  FURTHER  STUDY 

In  Chapter  III  it  was  noted  that  afterbody  shapes  other 
than  a  symmetric  foil  have  been  tried  without  significant 
improvement.   This  does  not  mean  that  this  approach  showed 
no  promise  for  possible  improvement  in  the  model,  but  only 
that  the  actual  shapes  tested  were  not  the  best.   It  is  felt 
that  this  model  can  be  improved  by  changing  the  shape  used 
for  the  jet  and  afterbody. 

As  reported  earlier,  the  crossflow  fluid  particles  near 
the  plate  apparently  see  a  cylindrical  shape  when  approaching 
the  jet.   As  they  pass  the  jet,  some  particles  are  entrained, 
some  join  the  turbulent  flow  within  the  wake,  and,  further 
out,  they  pass  as  if  the  jet  and  wake  together  are  a  solid 
object.   Just  downstream  of  the  jet,  along  the  crossflow 
axis,  there  is  an  area  of  high  entrainment .   This  causes  a 
high  local  pressure  drop  out  to  about  one  jet  diameter  for 
VR  =  8 .   Pressure  recovery  is  rapid  past  this  point.   This 
leads  to  the  idea  of  using  a  shape  such  as: 


This  shape  has  a  cylindrical  leading  edge,  the  maximum 
thickness  is  just  downstream  of  the  jet  center,  and  then  it 
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narrows  quickly  to  a  relatively  thin  tail.   The  total  length 
is  the  same  as  that  for  the  symmetric  foil. 

It  is  expected  that  this  shape  can  be  formed  using  the 
Joukowsky  transform  of  the  existing  model  by  varying  r  as  a 
function  of  the  angular  position  of  the  particular  point  at 
which  the  perturbation  velocities  are  being  calculated. 
Although  a  functional  relationship  which  will  work  over  a 
range  of  velocity  ratios  has  not  yet  been  f ound ,  it  is  pos- 
sible that  one  exists.   This  idea  might  also  be  pursued 
using  panel  (distributed  singularity)  methods. 

Another  area  of  uncertainty  during  the  course  of  this 
study  was  the  empirical  jet  trajectory  equation.   Reports  of 
experiments  to  determine  the  jet  axis  provide  a  wide  range 
of  possible  trajectories.   If  the  trajectory  equation  used 

in  this  model  is  accepted  as  correct,  the  trajectory  con- 

2 
stant  B  =  . 19VR   still  may  not  be  accurate  over  a  wide  range 

of  velocity  ratios.   In  particular,  it  is  suspected  that  it 

may  be  in  error  for  velocity  ratios  near  4  where  an  anomaly 

in  the  plate  pressure  distributions  occurs.   The  cause  of 

this  anomaly  is  not  clear,  but  it  is  suspected  that  this  is 

the  point  where  the  jet  momentum  is  so  low  compared  to  that 

of  the  crossflow  that  the  jet  reattaches  to  the  plate  within 

a  few  jet  diameters.   If  so,  the  jet  trajectory  would  be 

quite  different  than  that  which  is  currently  modeled. 

Finally,  the  entrainment  of  the  crossflow  fluid  into  the 

jet  is  implicity  accounted -for  in  the  formulation  of  the 
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horseshoe  vortices.   However,  the  pressure  drop  associated 
with  the  jet  entraining  the  surrounding  fluid  is  not  explic- 
itly determined.   This  model  predicts  a  zero  pressure  distri- 
bution at  very  high  velocity  ratios  (zero  ship  speed) . 
Further  improvement  of  this  model  in  accounting  for  entrain- 
ment  is  needed  if  it  is  to  be  used  to  predict  the  performance 
of  bow-thrusters  in  ships  with  no  way  on. 
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Figure  7.   Jet  axis  and  conforming  vortex  loop  system.  [7] 
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VR  =  4 


Figure  12. 


Plots  of  coefficients  of  pressure  along  the  free- 
stream  axis  used  to  determine  the  length  of  the 
foil.  Data  from  [5]. 
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experimental  findings.  [1] 
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APPENDIX  A 
DERIVATION  OF  THE  JOUKOWSKY  TRANSFORMATION 

2 
The  Joukowsky  transformation,  Z  =  x,    +  I   /z,,    is  used  to 

transform  a  foil  in  the  z-plane  to  a  cylinder  in  the  c-plane 

so  that  the  foil  can  be  treated  as  a  doublet  in  a  uniform 

flow.   The  terms  used  in  this  transformation  are: 

Z  -  the  complex  position  vector  of  a  point  in  the 

Z-plane  Z  =  x  +  iy 
?  -  The  complex  position  vector  of  a  point  in  the 

C-plane.   c  =  I   +  in 
l   -  A  positive  real  constant  equal  to  1/4  the  distance 

between  the  foci  of  the  foil, 
a  -  The  diameter  of  the  cylinder  in  the  5 ' -plane . 
m  -  The  distance  between  the  n  and  n'  axes.   This  is  a 

function  of  r. 
r  -  The  distance  from  the  forward  focus  to  the  leading 

edge  of  the  foil.   That  is  r  =  |c  -  2l\    where  -c  is 

the  abscissa  of  the  leading  edge. 
X,   -plane  -  The  linearly  transformed  c-plane  used  to  give 
the  foil  thickness,  c'  =  £*  +  in 

From  the  Joukowsky  transformation  and  Fig.  13  it  is  seen 
that: 


Z    =    .5(Z  ±  yj    Z2  -  Al2  )  (Al) 
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m  =  .5  (A  -  I)  (A2) 

z'    =    c  -(-m)  =  c  +  m  (A3) 

The  complex  potential  of  a  cylinder  in  a  uniform  stream 
is  given  by : 

biU')  =  U„(c'  +  a/c'J  (A4) 

The  complex  velocity  is  then  found  as  follows : 


du        doj        dc '       dc 

u»  +  u  -  1V  =  az  ■  ap    df     at 


do)        dco        d£*    /  dZ 
dZ  '     d?'      d?  /    ar 


then   substituting 


2 


2  2 


<|w   _   U     /  C    \         U         -a 


dz      -<<V     (?2  -  ,2> 


and   since    £      -    a   =    £    -    2, 


2 

dw    _    rT     /  5    \  (?  '    +   a) 


(A5) 


dZ         U~    t   ^ 
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Then 


rr       »      St 


(H 


U'    +   a) 


U    +    A) 


-    1 


(A7) 


FT       =   "^ 


/C\  U'    +    a)  I 

\v)     rr+rrj 


(A8) 
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APPENDIX  B 
COMPUTER  PROGRAM 


C  XXXXXXXXXXXX XXXXXXXXXXXXXXXX 

c  X  X 

C  X    THIS     PROGRAM    SIWULATCS     *     X 

C  X    FLUID    JET     INJ5=rrcr)    imtp       x 

C  X    A    UNIFORM    CROSSFLOW    AND       X 

C  X    FURTHER     DETCRMINES    THC  X 

C  X    PRESSURE    DISTRIBUTION    ON    X 

C  X    THE     SURFACE    FROM    WHICH  X 

C  X    THE    JET    IS     EJECTED  X 

C  X  X 

f  XX XX XX XX XXX XXXXXXXXXXXX XX XXX 

c 

C  CEFIMTICN    OF    TERMS 

c 

C      MR  AD    -  NUMBER  OF  RADIAL  POSITIONS  DN  SURFACE 

C  (INPUT) 

C 

C      MAMG    -  NUMBER  OF  ANGULAR  POSITIONS  ON  SURFACE 

C  (INPUT) 

C 

C  NPOINT    -    NUMBER    OF    ARBITRARY    POINTS    ON    SURFACE 

C  (MRAD*NANG) 

C 

C  UJ  -    JET    VELOCITY    (INPUT) 

C 

C  UM  -    MAIN    STREAM    VELOCITY    (IVO'IT) 

C 

C      O         ACTUAL  QtAMfTcR  hf  JET  nojeirF  (TMP'JT) 

C 

OR-    NON-OTMFNSIONA1     RADIAL    DOS!TTON    nF     SURFACF 

C  POINTS     (DIMENSIONED     AT    LEAST    NR AD )     (INPUT) 

r. 

C  THFTA       -    ANGULAR     POSIT]™,!     r>c    suPFACr     POTMTS 

r  (OIMFN'SIOMcD     AT    LEAST    MAMG) 

C 

C  XP  X-OOSTTION    OF    SURFACE    °0INTS 

C  (DIMENSIONED     AT     LEAST    NPOINT J 

C 

C  YP  Y-°nSITION    OF     SURFACE     POINTS 

C  (DINcNSinNED     AT     LEAST    NPOINT) 

C 

C      X         HORIZONTAL  COORDINATE  OF  J^T  AXIS 

C 

C  Z  VERTICAL    CCOPOIMATE   OF    JET       AXIS 

C 

C  ZINC         -    NON-OIMENSIONAL     INCREMENTAL     STED     SIZE    fpr    z 

C 

C      DELTAZ  -  Z-POSITION  OF  Ec F ECTI VE/VI R TUAL  SOURCE 

C 

C  HKSAV       -    CUMULATIVE     VCRTIf.  ITY     \LONG    T\-\<=    JCT 

C 

C  XI  NATUFAL    COORDINATE    TANGENT    TO    j  Fj  T 

C 

C  ZETA         -    NATURAL    COORDINATE    NORMAL    TO    J^T 

C 

C  ETA  -    NATURAL    COORDINATE    PERPENDICULAR    T0    XI 

C  AND    ZETA 

C 

C  ALPHA      -    ANGLE     BETWEEN    X-AXIS     AND    XI-AXIS 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


20 
30 


70 
C 


40 
50 

80 
C 


PAPEA  -    INCREMENTAL     PLATE     AREA 

TOTFCR  -    TOTAL    RESULTANT    FORCE    ON    PLATE 

X^AR  -    X-ROSITION    OF    RESULTANT    FQ-CP 

YT  -.NnRMAL!ZED    INDUCED    FORCE 

FLL  -    L/4    THE    DISTANCE     BETWEEN    TH^    FQCI 

OF    THE    FOIL     USED     FOR    THE    BLOCKAGE    EFFECT 

RAOJ  -    THE    DISTANT E    FROM    THE    FORWARD    FOCUS    OF    THE 

FOIL    TO    THF    LEADING    EDGE 
ZC°  -    THE    COMPLEX    VELOCITY    INDUCED     IN    the    UNIcORm 

FLOW    DUE    TO    TH=    PRESENCE    OF    THP    FOIL 


IMPL 
01  ^E 

X  VFL 
DIME 
CO  MP 

*Z?\  * 

DATA     IN 

READ 

REAH 

6C    FORM 

NPOI 

ZERO    0 

DO    3 

DO     2 

U(  I  , 

V{  I  , 

C°  (  I 

CO  NT 

CONT 

READ 

WR  I  T 

\^ELnCI 

VR     = 

RP  AQ 

Rc  AD 

FO  5M 

SET    CON 

31    = 

3HA 

34    = 

SET    PCI 

ELL 

RAD  J 

ALAM 

*4.  0* 
EM  = 
POSITI 
DCLT 
WR  IT 
WR  IT 
FORM 
FO  PM 

X' VEL 
FORM 

X »  D  T  A 
FORM 

ccfi  m 
pi   = 

TRAJ  FC 

3    = 


ICI 

NSI 

(45 

NSI 

LEX 

16, 

PUT 

(5, 

(5, 

AT( 

NT 

UT 

0    I 

0    J 

J) 

J) 

,J  ) 

IMU 

INU 

(5, 

E(  6 

TY 

(UJ 

(5, 

(5, 
AT  ( 

STA 

1.0 

=    0 

4. 

L     P 
=    6 


T    REALMS     (A-HtO-Z) 

ON    U(  50  ,S0)  ,V(50,=iO  )  ,W(  50,50)  ,R  {  50)  ,THETA{50  )  , 

O),CP(5  0,50) 

CM    XP(450)  ,Y°  (450) ,  Z°(450 ) 

ZED*16,!ST*16,ZCD* 16,ZSTP*16 , ZR*1 6, 
Z4*16tZHA*i6, Zw*16 


60) 
^0  ) 
14) 

=  N 
3T0 


=    0 

=    0 


NR  AD 
NANG 

RAD*NANG 
RAG-    MATRICES 
It  MR  AD 
1, NANG 

.0 
.0 
0.0 


40  ) 
,50) 
RATTO 
/'IM  ) 
70) 
70)  ( 
FLO. 
NTS 


JJ,UM,D 

UJ  ,UM,D 


ZINC 

7 (  I  )  ,I  =  L, NRAD) 

5) 

DSCD    IM    FCIL    3L0CKAGE 


CALCULATIONS 


.5 

\j 

AP 
.0 
.1 
.5 


ELL-** 

0.5* 

FN     PF 

AZ      = 

e  (  6 ,  a 

E(6,9 

AT(3F 

AT(  «0 

CC  I  TY 

ATI  »0 

yrrn^ 

AT     (' 

TION 

2.0* 

TORY 

0.19* 

AMETERS 
/VR 
2/f:LL 

*(2.*ELL  +  RADJ+DS0RT( ( 2.*ELL+RA0J  )**2- 
2)  ) 
(  AL  AM 
EFFE 
.6*DS 
0)  HF 
0  )  ZI 
3.4  ) 

•  ,  2  X  , 
-     '  , 

•  ,2X, 
S  A 130 
0'  ,2X 
OF  ^T 
ARS  IN 
CONST 
VR**? 


-    ELL  ) 
:CTIVE    SOURCE 
QRTCVR  ) 
LTAZ 
NC 

•JET    v/F.LOCTTY    =    •  ,    F8.  4,  2X,  •  UNI  F0D  m  •  , 

F8.4,2X,'J-T     DIAMETER     =      «.F3.4) 
•EFFECTIVE    SOURCE     IS,,F5.2,2X, 

IV^    TW    PL  "Tc  '  ) 
,'  Z    ST  EP    SIZE    =     ■ 

(1  .000  00000) 
ANT 


5.3  ) 
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c 
c 
c 


100 


c 
c 


CALOJLATICN    OF     INTERFERENCE    VELOCITY 

L    =    0 

CALCCN    =    r»FU0AT(MAN5  )    -    1.0 
INCREMENT    RADIUS     (INPUT    VALUES) 

DO     140     I    =    1  ,  NRAD 

RCALC    =    *(  I) 
INCREMENT    TH5TA    (5    DEGREE    INCREMENTS) 
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